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ABSTRACT 



The existing analytical solutions for gas lubricated stepped 
sector thrust bearings are reviewed. An experiment performed at high 
bearing numbers (123 to 165) using air at atmospheric pressure is 
described. Results are presented of tests conducted using a bearing 
optimized by the methods of Ausman and Wildmann for high speed operation. 
Most of the design optimizations are confirmed for high bearing numbers; 
the exception is load carrying capacity. The load capacity curve ob- 
tained takes the same form as that predicted by theory but is lower for 
any given bearing number by a constant percentage. This points to a 
need for re-evaluation of determination of load carrying capacity. 
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NOMENCLATURE 
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Bearing area (in. ) 

Radius ratio ■ r^/R 
Lubricating film thickness (in.) 

Film thickness in leading edge region (in.) 

Film thickness in trailing edge region (in.) 

Load parameter ■ W/p^ A 
Number of sectors 
Rotation rate (rpm) 

Pressure (psl) 

Compression ratio * h,/h. = ( A + h )/h„ = 1 + l/tx 

i 2 2 ^ 

Mass flow rate (Ibm./sec.-in.) 

Radial coordinate (in.) 

Bearing inner radius (in.) 

Bearing outer radius (in.) 

Time coordinate (sec.) 

Velocity of fluid (in. /sec.) 

Bearing load (Ibf.) 

Axial coordinate (in.) 

Minimum clearance ratio =* h 2 /A 
Step depth (in.) 

Angular width of groove between thrust pads 
(degrees or rad.) 

Bearing number “ — ^ — 

pH ^ 

Bearing number based on step depth = ' 

Modified bearing number = jfV 

Modified bearing number based on step depth = A-^^( I 
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y 

T 
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r 



Sector arc length (degrees or radians) 

Step ratio ** '^/y 

3 

Lubricant density (Ibm./ln. ) 

Coordinate angle (radians) 

2 

Viscosity (reyns = Ibf . -sec. /in. ) 

Arc length of leading edge region (degrees or radians) 
Angular velocity (rad. /sec.) 



Subscripts 

1 Refers to leading edge region 

2 Refers to trailing edge region 

a Refers to ambient condition 

r Radial 

© Circumferential 

A Refers to the step 

opt Optimized 
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1 . 



Introduction. 



In a rapidly expanding technology, which creates demands for machin- 
ery components that will operate in progressively more severe or unusual 
environments, it is soon found that "conventional" lubricants for moving 
or rotating parts will no longer serve. Because they possess several 
unique advantages, gases are a logical consideration as a lubricant when 
seeking solutions to many design problems. 

Gas lubricants are particularly attractive at extremely high or at 
very low temperatures. Since, in contrast to liquids, the viscosity of 
gases is relatively constant and even increases slightly with temperature, 
gas lubricants have an obvious application in very high temperature rotat- 
ing machinery. At the other end of the scale, elements such as nitrogen 
and helium retain their desirable lubricating qualities at cryogenic 
temperatures. Indeed, the "boiloff" from liquified gases is often used 
as the lubricant for expansion turbine gas liqulfiers. 

Another Important advantage of gas lubricated bearings is the extreme 
ly low power loss from friction at very high speeds. This, plus the fact 
that this power loss is relatively unvarying, dictates the use of gas 
bearings in precision applications such as high speed gyroscopes. The low 
friction loss at high speeds plus high temperature operation makes gas 
lubricated bearings the most desirable type for gas turbines and turbo- 
jet engines. 

In many areas gas lubricants can solve the problem of contamination. 
When pumping, compressing, or expanding high purity gases or refrigerants, 
the use of the process fluid as a lubricant will prevent mixing of the 
fluid with oils or greases. In food processing machinery, elimination of 
petroleum lubricants prevents contamination of foods. In nuclear reactors 
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it is desirable to eliminate oils, greases, and solid lubricants which 
break down or form radioactive isotopes when subjected to radiation. 

As with conventional bearings, gas lubricated bearings fall into two 
major classes, self-acting, or hydrodynamic bearings, and pressurized, 
or hydrostatic bearings. There exists a wide variety of both journal and 
thrust type bearings in each class. Self-acting bearings are, in general, 
the least complex, requiring neither separate gas supplies and piping nor 
fluid passages and orifices to complicate the design of a machine. 

Unfortunately, gas lubricated bearings have only recently emerged 
from the realm of laboratory curiosities. Currently most types are sub- 
ject to serious instabilities and require extremely close machining 
tolerances with very smooth finishes. Air film thicknesses are so small 
that thermal distortions are a major cause of failures. To gain a better 
understanding of the capabilities and problems of gas lubricated bearings, 
this project was undertaken. 

Because of an interest developed while the author was working at 
AiResearch Manufacturing Company on a summer industrial tour, this study 
will be confined to self-acting thrust bearings of the stepped sector 
type. There are several arguments which support this choice. 

Because the stepped sector thrust bearing (See Figure 1 for typical 
configuration) is a common, popular type noted for its simplicity and 
relative ease of manufacture, a large voluune of definitive theory and 
design data has been developed. The most recent and acceptable theory as 
yet lacks complete experimental verification. The recent development, 
at AiResearch, of a thrust bearing test rig capable of very high speeds 
affords an opportunity to obtain data in regions of operation previously 
unattainable. Arrangements were made for the construction of a second 
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Thrust Bearing Test Rig ML1606 at the Naval Postgraduate School in 
order to conduct this research. 
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2. Theoretical Considerations. 

In 1886 the English scientist, Osborne Reynolds, published the de- 
fining differential equation for thick film lubrication. Several years 
later, Kingsbury in the United States and Michell in Australia demonstrat- 
ed the applicability of the Reynolds Equation to thrust bearings by 
development of the tilting pad bearing. In 1918 Lord Rayleigh [1] found 
that for an infinitely wide, incompressible film, the optimum shape for 
a thrust bearing is the stepped configuration. However the step bearing 
can be optimized for only one expected operating speed or load but the 
tilting pad bearing automatically adjusts its pad angle to a new optimum 
as speed or load changes. Additionally, when side leakage is taken into 
account, the stepped bearing's load capacity markedly diminishes and its 
advantage over the tilting pad bearing disappears. For these reasons 
this bearing type was little used for liquid film lubrication. It was 
not until 1950 that Archibald [2] published the first analysis of a 
finite stepped thrust bearing using an Incompressible film. 

When gas is used as the lubricant, the step bearing becomes desir- 
able for two reasons. First, for a tilting pad bearing to retain its 
continuous optimum operation, the position of the pad pivot point must be 
continually changed. This is the effect of compressibility on the bear- 
ing design. A pivot is required that consists of a cam surface that is 
very difficult to machine. Secondly rvery close machining tolerances and 
finishes are required on gas bearings and since the tilting pad bearing 
is composed of a number of close-fitting parts, it is extremely costly 
to produce. In comparison the step bearing is relatively cheap. 

When the analysis takes into account compressible lubricants, the 
differential equation becomes nonlinear and several assumptions become 
necessary in order to obtain an analytic solution. Even though the 
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resulting solution is approximate It does establish a base for optimiza- 
tion of bearing geometry and permits correlation of theory to specific 
test bearings. 

Figure 1 depicts the geometry of a typical stepped sector thrust 
pad. The upper thrust plate (plane surface) rotates with constant angu- 
lar velocity while the lower (stepped) surface is stationary. Radial 
grooves of angular width 8 admit gas at ambient pressure to the lead- 
ing and trailing edges of each pad. The pad Itself, of sector arc 
length V , is divided into two sectors, one sector of length being 
recessed by an amount ^ relative to the other. The larger clearance 
In the recessed region Is designated h, , while the smaller, or mini- 
mum, clearance is designated h^ 

In cylindrical coordinates, the compressible lubrication equations 
for slider bearings take the form [3]: 
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where and are the radial and tangential components of lubri- 

cant velocity, is the radial mass flow rate per unit circumferential 
width and is the circumferential mass flow rate per unit radial width. 
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To obtain these equations, the following assumptions were made: 

(a) the flow is steady (does not vary with time); (b) the film thickness 



is at least two orders of magnitude smaller than r or r© ; (c) the 
axial, z, component of velocity and its derivatives are negligibly small; 
and (d) the gas lubricant behaves as a perfect gas and compresses and ex- 
pands isothermally within the bearing. 

By introducing the boundary conditions: 

V^(z=0) = 0 
Vp(z=h) = 0 
Ve(z=0) = r«b 
V^(z=h) = 0 

and integrating equations (1) and (2) the following expressions for the 
mass flow rates can be obtained: 



Q = - fJl JjL 
0 - Ip' 

L 2. a© 



(5) 

( 6 ) 



Ausman [4] introduced equation (4) and made the assumptions that (a) 

may be substituted for p , ignoring the effect of compressibil 

_2 

ity on the pressure gradient terms and (b) r in equation (6) may be re- 
-2 

placed with r^ where r, is some constant such that bR < r, < R where 
b is the ratio of inner to outer bearing radii = r^ /R. The first assump 
tion linearizes the resulting differential equation while the second makes 
the variables separable. Further, introduction of two dimensionless vari- 
ables, 

p = P- 



r =_£_ 



R 
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and substitution of equations (5) and (6) into equation (3a) results 
in the linearized, separable differential equation 



Jl jKf ll - ^ ^ 



(7) 



where 




4 R’*' 



( 8 ) 



Equation (8) defines the dimensionless parameter known as bearing number 
/V . The form of the solution is 



P=X(r)T(0) 

and results in an Infinite series solution which must be solved indepen- 
dently for each flow region* 

Boundary conditions state that the pressure equals the ambient pres 
sure at both inside and outside bearing radii, and that the pressure and 
flow rate are equal at the step when approached from either direction* 
These latter conditions are met by an iterative process which seeks a 
value for r^ which, when substituted into the two series solutions, re- 
sults in equal solutions at the step for both regions* 

The theoretical load carrying capacity for a bearing is found by 
integrating the expression for p over the bearing area* 

Ausman then optimized bearing geometry by independently varying 
step location, T/y , and number of pads, n, and solving for these para 
meters at maximum load capacity* He presented his results in the form 
of one table and one graph from which a designer can determine the com- 
plete optimum geometry for a given bearing size and expected operating 
point (bearing number. A- ) . In his optimization he assumed that the 
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optimum compression ratio, q, was 2.0 where * ^‘/ \r\ 2 _ ~ ( h^+ A )/h^ 

= 1 + ^ /ac . This, he found from trial solutions, was a good assump- 

tion over a wide range of bearing numbers. Note in Figure 3 (which is 
described later) that the solution for constant minimtim clearance is maxi- 
mum at a clearance ratio of 1.0 which corresponds to q = 2.0, 

Wildmann, e^. ^. , in a comprehensive study [5], obtained exten- 
sive experimental data at low bearing numbers which essentially con- 
firmed Ausman's theory in the regions the study undertook to investigate. 
Ausman's results were significantly simplified by the introduction of a 
modified bearing number 



A' = (10) 

This enabled him to eliminate the radius ratio, b, as an independent 
parameter and to present load capacity for optimum compression ratio 
(q = 2) as one curve plotted against bearing number. This curve is 
shown as Figure 2. 

As a result of experimental data, slight modifications were made to 
Ausman*s theoretical design parameters. New design parameters were pre- 
sented in graphical form so that design of a specific bearing could be 
undertaken without the necessity of extrapolating Ausman*s tabular re- 
sults. The first two design charts are similar to Figures 2 and 3 but 
reflect Wildmann*s correlation of theory. The remaining three are plots 
of optimum sector ratio, > and optimum sector arc length 

(degrees) and optimum ntimber of sectors, n, all plotted vs. modified 
bearing ntimber. The test bearing investigated in this study is designed 
using the parameters given by Wildmann. 
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One additional curve, computed from Ausman*s theory and presented 
in Wildmann's study, is shown as Figure 3 because of its usefulness in 
comparing experimental results with theory. This consists of a plot of 



linimum clearance ratio, ^ ;i / /S. 



Load Ratio = 



for a constant step depth, A . 



, vs . 

Load 

Optimized Load for q = 2 
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3. Description of Apparatus 

Trust Bearing Test Rig 

The primary research tool used in this experimental investigation 
is the AiResearch Thmst Bearing Test Rig, ML 1606, an extremely versatile 
device designed to measure the operating characteristics of gas lubricat* 
ed thrust bearings at very high rotation rates. Rotating speed and bear- 
ing load can be continuously varied and accurately measured while instru- 
mentation can be provided for measurement of air film thickness, torque, 
and temperatures at selected points. Design variations permit testing 
thrust bearings up to six inches in diameter and test bearings may be of 
any type or configuration. 

Photographs of the Thrust Bearing Test Rig and its associated instru- 
mentation are shown in Figures 4 and 5. A schematic of the arrange- 
ment of major components is presented in Figure 6. The rotating assembly 
consists of a hardened stainless steel shaft one inch in diameter with an 
integral 2% inch diameter thrust disc on the lower end and a small air 
driven turbine wheel on the upper end. This shaft is supported vertical- 
ly in two AiResearch ARPF compliant journal bearings. It is supported 
from beneath by the test bearing and separation of the back of the thrust 
plate from the upper body of the machine is maintained by an externally 
pressurized air thrust bearing. This configuration is capable of rota- 
tion rates up to 105,000 REM after irunning through the half -frequency 
whirl region encountered in the range between 20,000 and 50,000 RPM, 

The test bearing is supported beneath the thrust runner on the top 
face of a heavy steel piston which moves vertically and through a small 
rotation angle within a cylinder equipped with a hydrostatic air bearing. 

In the tests performed the test bearing was recessed into the face 
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of the piston and was glmballed by placing It on a 3/32 Inch thick sheet 
of gasket rubber. Restraint In torque was provided by securing a .020 
Inch X ^ Inch piece of piano wire In a groove on the piston face In such 
a way as to allow the end of the wire to extend Into the bottom of one 
of the air Inlet grooves on the test bearing. A photograph of the bear- 
ing secured In place Is shown as Figure 7. 

Air Supply and Control System 

A schematic of the air supply system Is shown In Figure 8. The 
system utilizes a 50 HP Westlnghouse 220 VAC, 3 phase, 60 hertz. Induc- 
tion Motor to drive a Chicago Pneumatic, Model P2CB Two Stage Air Compres- 
sor with Intercooler (max pressure = 450 psi). The output air Is cooled 
by a Hellflow shell and tube heat exchanger. Type 9AA14S. Water and oil 
are removed by a Wright -Austin Co., separator and float trap. These are 
followed by a Morse Diving Equipment Co., air filter and a standpipe con- 
taining a dessicant to further dry the air. Airflow to the receivers Is 
controlled by a Bailey Indicator Controller, Type K47, actuating a Bailey 
Dlaphram Valve Type % XSDS-9M. The compressed air is stored in two 

American Pipe & Steel Corp. , Vertical Air Receivers (each with a volume 
3 

of 117 ft . , max. pressure = 500 psi). Air is drawn from the top of 
both receivers and is distributed at outlets In a small manifold near the 
test setup. 

Figure 9 depicts schematically the air supply and flow control ar- 
rangement to the test rig itself. Figures 4 and 5 are photographs of 
the actual apparatus. Air is provided to the driving turbine through a 
one inch line fitted with a regulator and globe valve for rough speed con- 
trol. Fine control is achieved by a % in. line equipped with a regulator 
and needle valve. This line is placed in parallel around the rough control 
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line. A % in. line with a secondary filter and needle valve provides 
120 psig air to the hydrostatic thtrust bearing. Air is supplied to the 
hydrostatic piston journal bearing at 20 psig by a % in. line with regula 
tor and needle valve. A similar arrangement provides air to load the 
piston from beneath, however an atmospheric bleed is placed in this line 
to prevent pressure buildup due to low flow rates. Bourdon type pres- 
sure gages are located as shown in the schematic. Loading piston pres- 
sure is measured by a 14 inch mercury U-tube manometer. 

Electronic Instrumentation 

Figure 10 schematically describes the electronic Instrumentation 
used to measure rotational speed, air film thickness and temperature. 

Rotational speed of the test rig shaft is sensed by a Pace Engine- 
ering Co. SU-2 Pickup which is a very sensitive coil with a core of 
flexible iron straps. The coil, when placed in the vicinity of a rotat- 
ing shaft of ferromagnetic material, generates a small A.C. voltage of 
varying frequency as a result of the rotating magnetic field due to re- 
sidual permanent magnetism in the shaft. The Pace Pickup was placed as 
shown in Figure 4 and its output was displayed on the y-axis of a Hewlett 
Packard Model 130C Oscilloscope and a Hewlett-Packard Electronic Fre- 
quency Meter Model 500A. The output of a Donner Sine Wave Generator 
Model 1202 was displayed on the X-axis of the oscilloscope and its fre- 
quency was accurately measured by a Systron -Donner Digital Counter Model 
1013D. Since rotation rate is a predetermined figure for each data run, 
the rate it, preset on the sine wave generator and the speed of the rig 
is adjusted until a stationary elliptical lisajous figure is achieved 
on the cathode ray tube, indicating that the two frequencies are matched. 
This results in a speed measurement accurate to plus or minus one hertz 
or + 0.1% at 60,000 RPM. 
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Air film thickness measurements were made using Bently Proximitor 
distance detectors. The calibration and accuracy determination of this 
system is presented in Appendix A. During initial runs to check out the 
equipment it was found that heating and cooling of the distance probe 
(located in the center of the test bearing) was more extensive and more 
rapid than anticipated. This change in temperature causes a resistance 
change in the sensing coll and hence changes the output voltage for a 
given distance. Further, due to low friction, rundown time (time neces- 
sary to stop the shaft) was in excess of one minute and during this time 
the temperature and therefore the probe zero reading (thrust runner and 
bearing in contact), changed markedly. To correct this, two steps were 
taken. 

Two air nozzles made of % inch copper tubing were secured in the 
turbine outlet throat to direct air jets against the turbine blades op- 
posite to the direction of rotation. To actuate these, two spring load- 
ed, quick acting valves were installed. This system, seen in Figure 4 
stopped the shaft in an average time of 15 seconds. 

Two copper-constantan thermocouples were installed, one adjacent to 
the distance probe face and touching it, and the other in an ice-water 
bath as a reference junction. The output was recorded on a Mosely Auto- 
graf Strip Chart Recorder Model 7100B. The second channel of the re- 
corder was used as an event marker to record points at which data was 
taken. This arrangement is shown in Figure 10. 

These two corrective measures enabled the author to establish the 
probe temperature difference between the probe reading and the zero read- 
ing (thrust runner and bearing in contact) for each data point. The 
significance of this data will be shown in Appendix C - Experimental 
Uncertainty Intervals. 
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Test Bearings 



It was originally intended to test several thrust bearings of dif- 
ferent optimum configurations but time and funding limitations restricted 
the tests to one bearing. The bearing was chosen to correspond to the 
point in Ausman's Table I [4] at b = .20, VV = ® 160. The design 

parameters however, were taken from Wildmann [5]. These parameters are: 

b = .20 

Al= A.(l-t‘) = 153,7 
■? .135 

n s 3 

V = 118® 

= 16® 

The minimum film thickness was chosen to be 250. raicroinches so, for 
q = 2.0, the optimum step depth, A , is 250. microinches. The out- 
side radius, R, was chosen as 1.125 Inches. The final design based on 
these parameters is shown in Figure 11. Note the air inlet groove which 
is two degrees wide at the average radius. The material chosen for the 
bearing was "Graphmo" tool steel, a graphite -molybdenum steel noted for 
its wear resistance and low sliding friction. The bearing was manufactur- 
ed to specifications by B fie W Precision, Inc., Glendale, California. 

Inspection of the bearing by optical flat before running revealed 
a flatness of + 5.0 microinches across the bearing face. Profilometer 
readings by the manufacturer showed a surface roughness of 5 microinches 
(rms) . The step depth was measured on a Johansson Mikrokator Comparator 
and varied between 230 and 270 microinches with an average of 250. These 
readings were taken numerous times across all three steps. The surface 
of the recessed area was somewhat rougher than the raised portion but was 
not measured. 
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4. Experimental Procedure 

Data was obtained in such a manner that load vs. air film thickness 
could be plotted for a single constant speed. This was repeated for a 
number of different speeds to obtain a family of such curves. No data 
was taken in the region where half-frequency whirl was encountered to 
insure that the data was not affected by this instability. 

Test Rig Operation 

All runs were made with the hydrostatic thrust air set at 115 to 
120 pslg and the hydrostatic piston Journal air adjusted to 20 pslg. 

All starts and stops were performed with a test bearing load of 3.65 
pounds (3.0 inches Hg) . Prior to starting, the speed for the run (re- 
volutions per second) was preset on the digital counter by adjusting the 
sine wave generator to the same frequency (hertz) . The turbine was 
then started by quickly opening the rough control valve enough to bring 
the downstream pressure to about 5 psig. The bearing noise was monitor- 
ed continuously for whirl or rubbing by listening with an automotive type 
stethoscope. Rotational speed was adjusted while watching the llsajous 
figure on the oscilloscope and reading the frequency meter until the 
speed was slightly greater than the preset rate and steady. The tempera- 
ture recorder chart was then started and the bearing load slowly increas- 
ed, taking care that the temperature did not Increase at a rapid rate and 
cause excessive thermal distortion. As the bearing was loaded, the tur- 
bine air flow was increased to maintain speed. 

When the load for the particular data point was reached, the rig 
was allowed to run until the temperature reached a steady state or was 
increasing very slowly. Then the fine speed control flow valve was used 
to control turbine speed until the oscilloscope trace showed an exact 
frequency match. At this point, the distance probe output was read and 
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noted, the turbine air was secured, the braking jets were actuated, and 
the bearing load was reduced to 3.65 pounds. As the probe output was 
read, the event was marked on the strip chart recorder. At the instant 
that the shaft stopped and the thrust disc and bearing came in contact, 
the distance probe output was again read and the event marked. 

This procedure was repeated for each data point obtained, and each 
data point was repeated several times to insure correctness and to ob- 
tain a statistical average. The data obtained are presented in Tables 
1 through 4. 

Bearing Failure 

On the last run for data, a failure occurred and the thrust plate 
contacted the bearing face with subsequent rapid loss of air film and 
wiping of the bearing surfaces. The failure was initiated at a speed of 
72,000 RPM and a bearing load of 14.7 pounds. Just prior to failure a 
vibration noise began and increased in amplitude even though the load 
was reduced sharply. The center of both the bearing and thrust disc 
were wiped to a radius of approximately .35 inches which would seem to 
indicate failure due to thermal distortion. This is felt not to be the 
case, however, since neither temperature nor temperature rate were ex- 
cessive at the time, and the load, while high, was no higher than for 
several other runs when the air film thickness was even less. Rather 
it is felt that the speed used was a critical vibration frequency for 
some component, most likely the rubber sheet used to gimbal the test 
bearing. The appearance of the surfaces after failure merely indicates 
that slight thermal distortion was present at the time of failure. 

Data Reduction Calculations 

The data reduction calculations used to obtain the results presented 
in the next section are described in Appendix B. 
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5. 



Results 



The numerical results calculated from the data obtained are present- 
ed in Table 5. 

Figures 12 through 15 are plots of bearing load vs. air film thick- 
ness for each speed tested. All data points are plotted and extreme 
uncertainty limits are shown. Figure 16 is a composite plot of the same 
data using only the statistical averages for air film thickness. 

Figure 17 is a graph of dimensionless load parameter, L, vs. the 
modified bearing number, for a compression ratio, q »* 2.0. 

Perhaps the most useful form of presenting the results is seen in 
Figure 18. The ratio of actual load to optimized load at q » 2.0 is 
shown on the ordinate logarithmically while the clearance ratio, '=*< => 
is plotted on the abcissa. By the optimized load it is meant the opti- 
mized load from the theoretical curve in Figure 2 determined for the 
same bearing number as the particular data point under consideration. A 
clearance ratio of 1.0 corresponds to 'q = 2.0. All of the data points 
from Figure 17 are shown in Figure 19. 
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6. Comparison of Results with Theory 

In Figure 19, the results shown in Figure 17 are superimposed on a 
reproduction of Wildraann's Design Chart A [5], which predicts load carry- 
ing capacity at q - 2.0 for a given bearing number, J\. . Note that 

the results take roughly the same form as the predicted curve but the 
load capacity is signifigantly lower. 

In Figure 20, the curve from Figure 18 is replotted on a graph which 
shows both Ausman's theoretical prediction and Wildmann's experimental 
results. Note that the results of this study have not been normalized 
to a load ratio of 1.0 as was done in the study compiled by Wildmann. 

Two factors are immediately apparent. The first is that the slope 
of the curve obtained in the current study is essentially the same as 
Wildmann's correlation curve. Secondly, the new curve predicts a load 
capacity that is a constant fraction, 2/5 or 0.4 of that predicted by 
the design curves [f5] based on Ausman's work. 
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7 . 



Discussion of Results 



It is felt that the careful calibration of the probes used to meas- 
ure air film thickness, plus the temperature monitoring employed, re- 
sulted in reliable air film thickness measurements. Further, the repeti- 
tion of data points and use of average values reduced random error. This 
opinion is borne out by the smoothness of the load vs. film thickness 
curves and the consistency of the load ratio vs. clearance ratio curve. 

There are two sources of possible major error. The first is due to 
the near impossibility of constructing the perfect stepped sector thrust 
bearing upon which the theory is based. It was noted that the variation 
in the step depth was plus and minus 20 microinches or 12% of the speci- 
fied depth. This plus the rougher finish on the recessed area undoubted- 
ly affected the flow characteristics of the lubricant. The second seri- 
ous condition is thermal distortion. Even though attempts were made to 
allow temperature to reach a steady state, distortion was surely present 
as shown by the nature of the bearing failure which occurred. Thermal 
distortion would result in a nonuniform air film thickness; the film 
would increase in thickness as the radius at which it is measured is in- 
creased. It is not felt, however, that either of these conditions are 
responsible for the wide disparity between theory and experiment with 
regard to load capacity, particularly since the data is otherwise con- 
sistent. 
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8. Conclusions 



From the comparison of results with theory, it is concluded that: 

(a) Wildraann^s correlation of the variation of load ratio with clearance 
ratio is correct at high bearing numbers for the bearing tested; (b) 
Ausman*s theory correctly predicts the shape of the load vs, bearing 
number curve in the high bearing number region, and (c) the prediction 
of the magnitude of the load capacity is in error in this region. 

Ausman, in the preparation of his design table, performed separate 
optimizations for bearing geometry but extrapolated load capacity from 
curves computed for a central point in the table. He warned his readers 
that the accuracy of this extrapolation probably decreased with distance 
from the center of the table. The results of this test, of a bearing 
purposely chosen from an extreme comer of his table, confirm that Ausman *s 
intuition was correct. 



t 
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9. Recommendations for Further Study 

Now that it is established that there is a limit to the range of 
accuracy of Ausman's design parameters, additional bearings, optimized 
by the procedure used in this study, should be tested until the range 
of accuracy is established. 

A computer aided analysis of the load capacity equations should be 
performed over the full range of bearing numbers using the modified 
bearing number introduced by Wlldmann. The results should give some 
insight into the true nature of the relationship and could result in a 
revision of the optimized design parameters. 

Other types of thrust bearings, such as the spiral groove or her- 
ringbone patterns, await testing at high bearing numbers. The AlResearch 
Thrust Bearing Test Rig is ideal for this purpose. 
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TABLE 1 



Summary of 

RUN NO. 

14-1 

14-2 

14-3 

14-4 

14-5 

14- 6 

15- 1 

15-2 

15-3 

15-4 

15-5 

15- 6 

16- 1 

16-2 

16-3 

16-4 

16- 5 

16-6 

17- 1 

17-2 

17-3 

17-4 

17-5 

17-6 



Data -- 54,000 RPM 



LOAD 

(pounds) 


in.) 


6.59 


363 


6.59 


366 


6.59 


371 


6.59 


363 


6.59 


366 


6.59 


370 


8.62 


318 


8.62 


330 


8.62 


337 


8.62 


320 


8.62 


325 


8.62 


320 


10.65 


287 


10.65 


287 


10.65 


287 


10.65 


284 


10.65 


289 


10.65 


294 


12.67 


254 


12.67 


260 


12.67 


256 


12.67 


258 


12.67 


254 


12.67 


253 
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TABLE 1 (Cont) 



RUN NO. 


LOAD 






(pounds) 


( ^ in.) 


18-1 


13.69 


246 


18-2 


13.69 


246 


18-3 


13.69 


244 


18-4 


13.69 


251 


18-5 


13.69 


243 


18-6 


13.69 


250 
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TABLE 2 



RUN NO, 

1-1 

2-1 

2-2 

2-3 

2-4 

2- 5 
2-6 

3- 1 
3-2 
3-3 
3-4 
3-5 

3- 6 

4- 1 
4-2 
4-3 
4-4 
4-5 

4- 6 

5- 1 
5-2 
5-3 
5-4 
5-5 



Summary of Data -- 60,000 RHl 

LOAD 

(pounds) 

no data 

4.06 

4.06 

4.06 

4.06 

4.06 

4.06 

6.59 

6.59 

6.59 

6.59 

6.59 

6.59 

8.62 

8.62 

8.62 

8.62 

8.62 

8.62 

10.65 

10.65 

10.65 

10.65 

10.65 



h 



2 



in.) 



442 

452 

455 
466 

456 
455 
360 
363 
370 
380 
385 
394 
320 
313 
330 
333 
337 
328 

298 
294 
294 
296 

299 
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TABLE 2 (Cont) 





(pounds) 


2 

(t(in.) 


5-6 


10.65 


289 


6-1 


12.67 


256 


6-2 


12.67 


257 


6-3 


12.67 


256 


6-4 


12.67 


268 


6-5 


12.67 


261 


6-6 


12.67 


260 


7-1 


13.69 


243 


7-2 


13.69 


251 


7-3 


13.69 


253 


7-4 


no data 




7-5 


13.69 


253 


7-6 


13.69 


251 


7-7 


13.69 


253 


00 

1 


14.19 


240 


8-2 


14.19 


250 


8-3 


14.19 


241 


8-4 


14.19 


240 


8-5 


14.19 


243 


8-6 


14.19 


240 
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TABLE 3 



Summary of Data -- 66,420 RPM 



RUN NO. 


LOAD 


h- 




(pounds) 


2 

( in.) 


9-1 


6.59 


376 


9-2 


6.59 


380 


9-3 


6.59 


385 


9-4 


6.59 


385 


9-5 


6.59 


377 


9-6 


6.59 


374 


10-1 


8.62 


332 


10-2 


8.62 


347 


10-3 


8.62 


337 


10-4 


8.62 


339 


10-5 


8.62 


336 


10-6 


8.62 


337 


11-1 


10.65 


299 


11-2 


10.65 


305 


11-3 


10.65 


312 


11-4 


10.65 


295 


11-5 


10.65 


295 


11-6 


10.65 


298 


12-1 


12.67 


263 


12-2 


12.67 


274 


12-3 


12.67 


264 


12-4 


12.67 


274 


12-5 


12.67 


268 
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TABLE 3 (Cont) 



RUN NO. 


LOAD 


^9 




(pounds) 


2 

(^in.) 


13-1 


14.70 


250 


13-2 


14.70 


247 


13-3 


14.70 


240 


13-4 


14.70 


247 


13-5 


14.70 


243 


13-6 


14.70 


246 
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TABLE 4 



Summary of Data — 72,000 RPM 



RUN NO. 


LOAD 

(pounds) 


^2 

( in.) 


19-1 


8.62 


337 


19-2 


8.62 


353 


19-3 


8.62 


346 


20-1 


10.65 


301 


20-2 


10.65 


301 


20-3 


10.65 


301 


21-1 


12.67 


265 


21-2 


12.67 


267 


21-3 


12.67 


277 


22-1 


13.69 


258 


22-2 


13.69 


260 


23-1 


14.70 


253 


23-2 


14.70 Bearing 


failure. 
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Tabulation of Results 
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TABLE 5 (Cont) 

Run No. Speed Load Avg. Bearing No, L L 

(rpm) (Ibf.) (^in.) (A^) 
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SECTION A-A 
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NOT SHOWN TO SCALE — CLEASMCtS GREATLY EXMOCBtAIED 



FIGURE 1* TYPICAL STEP BEARING GEOWTRY 
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0 15 30 45 60 75 90 105 120 135 150 165 180 

Bearing Nuaber CA^ 

FIGURE 2. THEORETICAL LOAD CAPACITY FOR COMPRESSION RA 10 = 2,0 



0ptimi2ed Load for q‘^2 



u 

o 

•o 




FIGURE 3. 



•2 .4 .6 .8 1.0 1.2 1.4 1.6 1.8 2.0 

Clearance Ratio = 

THEORETICAL VARIATION IN LOAD RATIO AS A FUNCTION OF 
CLEARANCE RATIO DETERMINED FOR A « 30 
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FIGURE 4„ AIRESEARCH THRUST R^\R1^G !LST RIG, ML1606 



A7 



OSCILLOSCOPE 




^8 



FIGURE 5. OVERALL VIEW OF TEST APPARATUS 




V 




FIGURE 6. MAJOR COMPONENTS OF THRUST BEARING TEST RIG 
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PliiTOn FAf' 




K,rt4K 7. IKSl Bi.ARlNG IN PLACE ON PISTON rACT 




MR 

ZOOUN<i WATER 




FIGURE 8. SCHEMATIC OF AIR SUPPLY SYSTEM 




52 




STRIP CHART 

recorder 
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FIGURE 10. ELECTRONIC INSTRUMENTATION SCHEM.ATIC 




FIGURE II. TEST BEARING DESIGN 



Bearing Load, W ( Ibf 




Air Film Thickness, in.) 



FIGURE 12. BEARING LOAD vs. AIR FILM THICKNESS - 54,000 RPM 
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220 240 260 280 300 320 340 360 380 400 420 440 460 

Air Film Thickness* (^in.) 



FIGURE 13. BEARING LOAD vs. AIR FILM THICKNESS - 60,000 RPM 
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Bearing Load, W ( lb£ 




220 240 260 280 300 320 340 360 380 400 420 440 460 



Air Film Thickness, (^in.) 



FIGURE 14. BEARING LOAD vs. AIR FILM THICKNESS - 66,420 RPM 
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220 240 260 280 300 320 340 360 380 400 420 440 460 

Air Film Thickness, (^in.) 



FIGURE 15. BEARING LOAD vs. AIR FILM THICKNESS - 72,000 RPM 
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220 240 260 280 300 320 340 360 380 400 420 440 460 

Air Film Thickness, (^in.) 



FIGURE 16. BEARING U)AD vs. AIR FILM THICKNESS AT VARIOUS SPEEDS 
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Load Parameter 




no 120 130 140 150 160 170 180 190 

Bearing Number (-/O 



FIGURE 17. LOAD PARAMETER vs. MODIFIED BEARING NUMBER FOR q = 2 
(TEST BEARING) 
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Clearance Ratio = 

FIGURE 18. LOAD RATIO vs. CLEARANCE RATIO (TEST BEARING) 
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15 30 45 60 75 90 105 120 135 150 165 180 

Bearing Number (A^ 

FIGURE 19. THEORETICAL AND EXPERIMENTAL LOAD CAPACITY 




Clearance Ratio = 

FIGURE 20, THEORETICAL AND EXPERIMENTAL LOAD RATIO vs. CLEARANCE RATIO 
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APPENDIX A 



DISTANCE PROBE CALIBRATION 
Selection of Measuring Systems 

The thickness of the air film between the test bearing and its 
thrust runner was expected to vary between 150 to 400 microinches. It 
was felt that significant data could be obtained with a measuring techni- 
que which could discriminate within plus or minus 10 microinches during 
actual operation. Based on descriptions of apparatus found in the liter- 
ature, two types of sensors should be considered a capacitance probe 

which would use the air film to be measured as the dielectric between 
two parallel plates, and the magnetic, or reluctance, type which senses 
the proximity of a conductor by the extent to which it disturbs a mag- 
netic field surrounding a sensing coil. 

Although it is capable of high accuracy, the capacitance type was 
eliminated because it depends on the dielectric constant of the air film. 
It was felt that changes in temperature and moisture content would un- 
predictably alter the dielectric properties of the air film. 

Bently Proximitors manufactured by the Bent ly-Nevada Corp.were chosen 
because of their high potential accuracy, convenient mounting arrange- 
ments, and low cost. Since they are of the magnetic reluctance type, 
their measurements depend upon the nature c ' the conducting surface being 
sensed. Variations in the properties of nonconducting materials be- 
tween the probe and the measured surface have no measurable effect on 
the accuracy obtained. Temperature of the probe itself was found to 
have an effect, however, even though silver-wound coils, advertised as 
the least temperature sensitive, were used. 
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Enhancement of Accuracy 

The selected probes were Intended to measure a range of 0 to .010 
Inches and to resolve to + .0001 inches. To attain the necessary ac- 
curacy (+ 10 microinches), the resolution had to be improved tenfold. 

According to the manufacturer, the accuracy and resolution is direct- 
ly proportional to the stability of the input voltage. Accordingly, an 
ultra-stable power supply meeting the voltage and current requirements 
was purchased from Systems Research Corporation. Measured ripple and 
noise of their Model 3566 Power Supply was found to be 50 to 80 micro- 
volts peak to peak and the output varied no more than 1.0 millivolts in 
a 24 hour test. 

To read the output voltage accurately, a digital voltmeter. Electro 
Instruments Co. Model DVX-500, capable of measuring to the nearest 0.1 
millivolts was obtained. 

Probe Calibration Methods 

Photographs of the probe calibration rig used are shown in Figures 
A1 and A2. The apparatus consists of a Johannson Mikrokator Model S 
510E-7 Comparator altered for use as a calibrator. The Mikrokator dial 
indicator seen in Figure A2 is accurate to the nearest one millionth of 
an inch and is thus capable of establishing the actual resolution of the 
distance probes. The comparator was altered by providing a rigid clamp 
fixture to hold a distance probe over the movable table, and by making 
arrangements to adjust the position of the dial indicator head by means 
of a fine adjustment mechanism salvaged from a microscope. 

The surface to which the probes were calibrated was the surface of 
a disc made indentlcal to the thrust runner actually used in the thrust 
bearing test rig. This was done by cutting the disc from the same billet 
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of steel and having it follow the shaft with thrust runner through the 
same machining, heat treating, grinding and lapping processes. 

Provisions were made for checking the effects of temperature on the 
distance measurements. The ambient temperature was accurately recorded 
before, during, and after each calibration. The calibration rig was al- 
lowed to "soak" for at least three hours at any given temperature to elim- 
inate any transient expansion or contraction. 

A temperature sensitivity test was performed to insure that tempera- 
ture changes in readout were not due to mere distance changes as a result 
of expansion or contraction of clamping fixtures of the calibration rig 
itself. A probe was clamped in a vertical position with its sensing sur- 
face up. The test surface disc was allowed to rest above the probe i:ace 
with a separation of .015 inches maintained by a thin sheet of Pyrex 
glass. The coefficient of themal expansion, , of Pyrex glass is 
2.0x10 ^ ® temperature change of lOOF, for instance, the 

glass thickness ,a , will change by 

A ^ ^ J 

= (2.0xl0‘^)(.015)(100) 

=3.0 microinches 

This arrangement was then subjected to heated air blown from a portable 
heat gun while the output reading was watched for deviation. 

Results 

In the temperature sensitivity test, raising the temperature ap- 
proximately 50F caused a change in reading of about 100 millivolts. This 
represents a relative distance change of about 140 microinches and con- 
clusively confirmed suspicions that the probes themselves are affected by 
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temperature. 



Calibration curve data are presented in Tables A1-A3. The resulting 
calibration curves are presented in Figure A3. Note that data was taken 
only in a small linear region of the probe's effective range since larger 
displacement measurements were not anticipated. The slope of the curves 
was the same in every case: 



Slope » 1.4118 



microinches 

millivolt 



Discussion 

From Figure A3 it is immediately apparent that although temperature 
changes affect the position of the curve, only the zero, or reference, 
point is affected. The slope of the curve remains constant over the ex- 
pected range of measurements. Accordingly, the experimental procedure 
during actual runs Included the determination of a zero point (bearing 
and thrust runner in contact) immediately after each run while the rig 
and probe were still at operating temperature. Measurements were then 
taken from this zero point using the constant slope of the calibration 
curves. 

To determine data scatter and accuracy, the data was originally 
plotted on 20x24 inch graph paper with 10x10 squares to the half-inch. 

The mean deviation from the straight line plots was 0.96 microinches. 

Thus by accuracy and stability enhancements, together with careful cali- 
bration, the resulting data may be considered reliable within + 1.92 micro- 
inches, or twice the mean deviation. This is better than five times the 
necessary accuracy and fifty times better than the manufacturer's intended 
accuracy. 
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TABLE Al 



Run No, 



DISTANCE PROBE CALIBRATION DATA 
68,6°F 

Relative Distance 
(in.) 

400 

379 

337 

296 

268 

204 

170 

121 

66 

3 

392 

359 

320 

271 

216 

168 

111 

63 

9 



Probe Output 
(mv) 

6069 

6035 

6000 

5968 

5945 

5900 

5875 

5838 

5800 

5756 

6056 

6031 

6005 

5968 

5933 

5896 

5857 

5820 

5777 
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TABLE A2 



DISTANCE PROBE CALIBRATION DATA 



75.4“F 

Run No. Relative Distance Probe Output 

(in.) (rav) 

3 398 5943 

350 5907 

302 5874 

251 5836 

193 5795 

139 5759 

95 5724 

52 5695 

4 398 5941 

344 5898 

310 5872 

229 5815 

176 5777 

114 5731 

62 5689 

0 5648 
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TABLE A3 



DISTANCE PROBE CALIBRATION DATA 



49.4°F 

Run No. Relative Distance Probe Output 

(in.) (rav) 

5 396 6242 

366 6219 

307 6175 

254 6134 

234 6123 

177 6082 

217 6037 

68 6005 

0 5956 

6 395 6207 

360 6186 

310 6154 

251 6117 

207 6088 

158 6054 

123 6028 

87 6004 

51 5981 

18 5956 

0 5942 
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FIGURE Al. DISTANCE PROBE CALIBRATION APPARATUS 
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FIGURE A3^ DISTANCE PROBE CALIBRATION CURVES 



APPENDIX B 



DATA REDUCTION CALCULATIONS 

The results of measurements shown in the Summary of Data (Tables 
1 through 4) were calculated from the raw data. Examples of these calcu- 
lations follow. 

Since the speed measurement is indicated by a frequency (hertz) and 
1 hertz = 1 revolution per second 

speed in revolutions per minute is computed by multiplying this value by 
60. For example: 

N = 900 hz. X 60 = 54,000 rpm 

The bearing load is indicated by the height of a column of mercury 
in inches. The manometer reading is converted to psi and multiplied by 
the piston area. The weight of the piston and test bearing assembly is 
subtracted from the result to obtain the actual bearing load: 

2.036 in. Hg = 1.0 psi 
Piston Area = 10.317 in. 

Weight of Piston Assembly = 11.65 lbs. 

Therefore, for a reading of 4,0 in. Hg, 



W = 



4.0 X 10.317 
2.036 



11.65 = 8.62 lbs. 



Air film thickness is determined from the difference (in millivolts) 
between the voltmeter reading at a particular speed and load and the 
reading immediately thereafter taken with the bearing and thrust disc in 
contact. This difference is multiplied by the constant slope of the 
probe calibration curve (Appendix A) to obtain air film thickness. From 
Run No. 14-1, 
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h = (5864 - 5607)(1.4118) 



= 257 X 1.4118 = 363 microinches 

Subsequent analysis of the results is performed using the dimension- 
less parameters, > snd Load Ratio (L/L ). The calcu- 

lation of these parameters is described in the following paragraphs. 

The bearing number, \ ^ > is 



where constant values are 

R 
b 

and 

P 

For Run No. 6, 

P 

and 

a' = (6) (28 X 10'^°) (628 3 ) (1.125^) ^ 

(14.89) (25 X ^ 

= 137.8 

The load parameter, L, is given by 




where 

A = (R^ - rj)rrx (354/360) 

= (1.125^ - .225^) 7- (354/360) 
= 3.755 in. 



28 X 10 reyns 



= 1.125 in. 

= 25.0 X lO’^ in. 

= 0.2 



_ / 






! • X 



I / 



> £,t 



= barometric pressure (in. Hg) x .4912 ^ 



'J- 



1000 (27t) = 6283 rad. /sec. 
30.31 X .4912 = 14.89 psia 
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The fraction 354/360 accounts for the fact that the air inlet grooves 
are not part of the load carrying area. Again for Run No. 6, 



L = 



12.67 



(14.89)(3.755) 

= 0.227 

The clearance ratio, c< = , is computed from the average 

film thickness for a given run. Using Run No. 6, 



and 



h 2 (avg) = 259.7 microinches 



^ = 1.04 



250 



The optimized load for q = 2 is found by consulting Figure 3 for 
the particular bearing number, • For example if 

137.8 
L = .628 

and 

Load Ratio = 



.227 



.628 



= .361 
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APPENDIX C 



EXPERIMENTAL UNCERTAINTY INTERVALS 



Reference ^9' gives the uncertainty interval for a result as 



/ 




v/aero w , w w are the uncertainty intervals of the variables 

making up the result, R is the expression for the result, and v , ^ > 

v,^ are the variables making up the result. 

The expression for bearing load is 

A, - 

where 

p = load pressure (in. Hg) 

A^ = piston area (in.) 

= weight of piston & test bearing assembly (Ibm.) 

The uncertainty interval for bearing load is 



Dividing through by W to nondimensionalize the equation results in 





For 



Pt “ it (20 to 1) 



A^ = 10.317 + .001 in. (20 to 1) 




.01544 
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and 



w 



V7 



= (10. 65)(. 01544) = + 0.164 Ibf. 



Air film thickness is computed in the following way; 



hj . (Vj - V^)S 



where 



= probe output voltage at speed and load (mv) 
V 2 = probe output voltage at zero point (mv) 

S = slope of probe calibration curve (A| in,/mv) 



One additional factor must be taken into account. Inspection of the probe 
calibration curves shows an average temperature sensitivity of 14 mv/°F 
at a constant distance between probe and measured surface. Examination 
of all temperature data taken shows that in only one out of ten readings 
did the probe temperature differ at speed and load from the zero point 
by more than 1,25 °F. An actual temperature trace from the strip chart 
recorder is shown as Figure C~l. Defining an additional error, 



It is felt by the writer that repetition of each data point an average of 
six times halved this uncertainty interval at the very least. The inter- 
val used is 



w^^ = uncertainty interval due to a difference in prcbe 

temperatures 



= (14.0)(1.25) 



= + 17,5 mv 



w 



=+8,75 m V 



(10 to 1) 



Therefore for Run No. 5-2 using the previous method, 




where 
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= 5788 + 2 mv 



(20 to 1) 



= 5580 + 2 mv (20 to 1) 



= 1.4118 4- .005(iin./mv (20 to 1) 



and 



• -r— V4- I L'' 

: /'I ~ 1.1 -.rc/oj V^o; 

Wi^ = 294 (.0422) 

=112.41 microinches 
The expression for bearing number is 

« c '\'.o R' 



.V" 



^ ■ £1 
OS 



‘ u 



f \ /V' 






and by similar means 









A.' 






(^j'^ (^) - 2 (. 



where for Run No. 5-2, 



/( = 28x10 reyns 

cO = 6283 + 6.283 rad. /sec. 
p, = 30.31 + .01 in. Hg 
L. = 250 + 10 microinches 



and 



. 

'' i.'Wd) 



i ^ 

A 



(20 to 1) 
(20 to 1) 
(20 to 1) 
(20 to 1) 



/37.8 



/ L\ 4 



■ /i-2!iVa 



OI 



31 J 



:)■- ^ (^r. 



w. = 137. 8(. 0669) 

= + 9.22 

Similarly for Run No, 5 (averaged points) 

/ 

/a 
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4 - 



Wl. 






W'l 










Woc 

TTF 



/iHLV 

[ I 



(-T 

\%S£> ) _ 



/z 



= 1.18 (. 0581 ) 

= + .069 

The foregoing are examples of the determination of uncertainty 
intervals for each data point. These intervals are represented on all 
graphical presentations of data as short lines through the data points in 
the following manner: 
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FIGURE Cl 



PROBE TEMPERATURE AS A FUNCTION OF TIME 
FOR RUN NO. 17-4 
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